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SIMPLIFIED METHODS FOR INTERPRETING THE EFFECT 


OF TRANSFER-FUNCTION ZEROS ON THE 
TRANSIENT RESPONSE OF AIRCRAFT 


Reiner Onken* 
Ames Research Center 


SUMMARY 


Two simple methods are outlined for evaluating the effect of transfer- 
function zeros on the system time response. The pole effects can also he 
evaluated. These methods are useful for simplified analysis or creating 
design criteria in terms of desirable regions of pole-zero locations. 

The type of transfer function studied is limited to those of linear sys- 
tems. Corresponding to ordinary longitudinal or lateral aircraft transfer 
functions , the denominator polynomial is of fourth order and the numerator of 
third order at most. 

With the longitudinal motion of the aircraft- as an example s the methods 
are used in the evaluation of optimal regulator control with respect to a 
particular performance index structure. 


INTRODUCTION 


The dynamic behavior of a flight vehicle usually can be described by a 
system of linearized differential equations with constant coefficients or by 
the corresponding set of transfer functions. The use of feedback control in 
multivariable systems makes it possible to modify not only the denominator of 
the transfer functions, but also the numerator. This property is desirable 
for it leads to the additional possibility of improving the vehicle's dynamic 
behavior. As an extreme case, nonminimum phase effects (numerator roots with 
positive real component), which occur rather often in aircraft transfer func- 
tions, could be ruled out or at least be moderated by appropriate feedback 
control. 

Therefore, some work was directed toward the development of methods for 
the investigation of the numerator roots based on illustrative frequency 
domain techniques (refs. 1-U). These techniques provide an estimate of 
changes in the numerator and denominator roots of the transfer functions as a 
consequence of feedback loop gains. To establish design criteria in terms of 
appropriate root locations , the effect of the zero locations on the dynamic 
behavior, that is, t ime response, should be known. The time response appears 
*National Research Council Associate. 



to be most suitable for evaluations of system behavior, since the necessity of 
an analytical formulation of the performance index can be avoided. System 
time response may be obtained by resolving the transfer functions into partial 
fraction series and taking the inverse term by term (residue calculation). 

This procedure is well known, but it does not offer a means of showing trans- 
parently what changes in the time response can be expected from certain pole- 
zero modifications. Especially for multivariable systems, a more effective 
procedure is required. Various pole-zero configurations have to be evaluated 
to determine the effect of modifications on the corresponding time histories. 

This report discusses simple methods that show the effect of zeros and 
poles on the transient response so that root locations can be assigned for 
certain desirable transient response behavior. Two basic methods will be 
described: The first is based on the assignment of phase and magnitude of the 

residues to the root locations; the second, which is easily applicable to air- 
craft control systems, is based on a specific expansion of the transfer func- 
tions, where each expansion term itself has a known transient response. 

The application of both methods is restricted to a controlled element 
transfer function with a fourth-order denominator and at most third-order' 
numerator. This typifies either the longitudinal or lateral directional 
equations of motion of an aircraft. 

Implementation of the methods is discussed for the evaluation of an air- 
craft optimal control system. The optimal longitudinal motion regulator with 
respect to a particular quadratic performance index is evaluated for differ- 
ent relative weighting between the states and the controls. 


RESIDUE INTERPRETATION METHOD 


Transfer Function and Transient Response Residues 

The transfer functions investigated here are of the fourth-order denom- 
inator and at most third-order numerator type. This is the structure of the 
main open-loop aircraft longitudinal motion transfer functions as well as the, 
corresponding closed-loop transfer functions with constant feedback gains. If 
u c is generally taken as the input and x as the output, the transfer 
function becomes 


F(s) 


x(s ) 
u c (s) 


K 


( — P 1 ) (“P2 ) ( ~P 3 ) (“P4 ) 
(s-pi ) (s-p 2 ) (s-p 3 ) (s-p 4 ) 


(s— Z! )(s— z 2 ) (s— z 3 ) 
(-z 1 )(-z 2 )(-z 3 ) 


( 1 ) 


Disgarding the case of multiple roots, in the frequency domain the response to 
a step input can be written as a partial fraction expansion: 
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The inverse is 
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The partial fraction expansion coefficients , weighting the different modes 
involved, are 


Thus 
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Normalizing the transient response with respect to the steady state (K = l) 
gives C Q = 1. The remaining residues become 


where 
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The entire residue expression has been separated into two main factors C^p 
and Ci z , similar to each other, where one part (Cip) consists only of pole 
terms, while the other part (Ci z ) contains all zero-affected terms. Thus, Ci z 
covers the whole effect of the zeros on the residue. 


The zero-affected portion Ci z consists of products of fractions whose 
factors are complex numbers in terms of the poles and zeros . They could be 
drawn in the root plane as vectors from the different zeros to the pole p^ 
and the origin. For two zeros, four vectors would have to be investigated for 
each residue. This is not convenient, and a different representation of the 
zero effect on the residues is desirable. Furthermore, the representation 
should illuminate zero regions of similar effect. Such a representation is 
developed in the following sections. 

The pole-dependent portion C^p has a reciprocal structure to that of 
C^ z . Thus, the influence of the poles on the residues can be considered in a 
way similar to that of the zeros . 
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Once the effect of the zeros and poles on the residues is known, the 
time vector presentation is a helpful tool for putting together the various 
portions of equation (3) to determine the overall transient response (refs, k 
and 5 ). There are two oscillatory eigenmodes for aircraft longitudinal con- 
trol: the low-frequency phugoid mode and the fast short-period mode. Conse- 



Figure 1.- Residues and time vectors 
for given pole combination. 


quently, there are two time vectors 
5q and X3 at the time t = 0 (see 
fig. l) in phase with Cj and C3 and 
with a magnitude of twice Cj and C3. 1 
They rotate with time corresponding to 
their eigenmode frequency, while 
decreasing in magnitude as a result of 
their eigenmode damping. The tran- 
sient response for any time t , due to 
the residues C* through C4 in equa- 
tion (3), can he read from figure 1, 
as the sum of the real portions of the 
time vectors. Thus, the purpose here 
is to determine the change of the 
residues (time vectors) in phase and 
magnitude introduced hy the zero- 
dependent numerators of the transfer 
functions . 


Phase Loci of the Residues 

As already indicated, the residues must he considered in general as con- 
jugate complex quantities, and they must .he described hy their phase and mag- 
nitude for later use with the time vector representation. Therefore, the link 
between the zero location and the residues is required before the overall cor- 
respondence between the zero locations and certain desirable transient 
response behavior can be established. Thus , the phase contribution to the 
residues due to the zero locations is presented here in terms of phase loci in 
the root plane for the following cases : 

1. Pair of conjugate complex zeros 

2. Pair of real zeros 

3. Three zeros 

Pair of conjugate complex zeros — The portion of the residue that is 

influenced by the zeros is known from equation (6) as 


C. 

1Z 




( 7 ) 


1 The residues C2 and C4 are implicitly included in Xj and X 3 . They do 
not appear explicitly, because they are given as the conjugate complex numbers 
Ci and C3 . 






Figure 2.- Definition of *P\ and ^2 • 


The phase angle 4>i represents 
the change in phase due to the zeros . 
From the four vectors involved, as 
shown in figure 2, the phase angle 
becomes 

= *1 + *2 ~ = *1 + *2 ( 8 ) 

with ^ = 0 for pairs of zeros. To 

determine the loci of zeros that have 
equal phase changes , consider 


tan <f> . 
1 


const 


tan Vi + tan ^2 
1 - tan • tan <^ 2 


(9) 


If the definition of and <P 2 is substituted from figure 2 into equation (9), 
it can be shown that those loci are circles of the form 


where 
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These loci are illustrated in the 
upper half plane in figure 3, where it 
becomes evident that all circles 
intersect the pole p^ and have their 
center on the real axis . They degen- 
erate to a vertical straight line 
through pj_ if $i is either 0, 2 n , 

. . . , or ir , 3tt , . . . If the por- 
tions of the circles on the left side 
of pi correspond to the phase angle 
<f>i , the right-hand portions corre- 
spond to 4»i + ir. For each circle, 
the angle 4^ can easily be 
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determined as the angle between the real axis and the line connecting the 
center of the circle with p^. 

Thus, each pair of conjugate complex zeros lying on one of the circles 
causes the same phase shift on the residue Cj. as any other nair of 

zeros on the same circle. The intersection points of the circles with the 
real axis correspond to double zeros. This yields the link between the case 
of conjugate complex zeros and that one of a pair of real zeros. 

Pair of real zeros — -The phase angle f° r the ca se of real zeros is 

also given by equation (8) (see fig. k ) . The loci of pairs of real zeros sat- 
isfying the condition of equal phase changes are 


where 


( °z - 



Drawing it in a a z ,c o z plane (the 
lower half plane in fig. 3) , we 
obtain a set of rectangular hyperbolas 
with their asymptotes inclined by 45° 
to the real axis. For corresponding 
values of 4>i , the hyperbolas inter- 
sect the real axis at the same points 
as the circles for the case of conju- 
Figure 4.- Definition of 'fi p and gate complex pairs of zeros. Thus, 

for pair of real zeros. considering both the base of conjugate 

complex zeros and that of a pair of 
real zeros , a certain range in the 

residue phase angle will correspond to an area in the root plane, bounded by 
the appropriate loci curves . 

Three zeros — If the number of zeros is extended from two to three, the 
equation for the residue portion affected by the zeros Ci z is slightly 
changed to 
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Figure 5.- Phase shift A<f>£ 
to third zero Z3. 


due 


Considering the phase loci, there is 
the same family of curves in the root 
plane as already derived, with an 
additional phase shift A<J>£ due to 
the third zero ; that is , 


4 > i = 4 > i + A<f>^ (13) 

1 (two zeros) x 


where A^ can he easily derived for 
the third zero moving along the real 
axis , as shown in figure 5 • There is 
a discontinuity in A 4^ of tt, when 
Z3 passes the imaginary axis. 


Magnitude Loci of the Residues 

As outlined previously, a mapping of the magnitude, as well as the phase, 
of in the root plane by means of locus curves is required.. Following the 

same line of reasoning as for the. phase effect of the zeros on one 

obtains locus curves of higher order, which are too complicated for simple 
geometrical interpretation (see appendix A) . Therefore, the expression for 
|Cj[ Z | will be split so that it can be discussed partially by simple geometri- 
cal means. Again, the cases of conjugate complex pairs of zeros, two real 
zeros, and three zeros are discussed. 


Pair of conjugate complex zeros — Following equation (7), the magnitude of 
C^ z is split into 


where 
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|C iz l ■ c ! * c 2 
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Therefore , the two portions Ci and c 2 are discussed separately and then their 
product |Cj_ z |. The loci curves in the root plane for both Ci and c 2 again 
are circles (fig. 6) given, respectively, by 
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Figure 6.- Magnitude loci of c \ . 
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Both families of circles are very similar: the circles for C 2 represent the 

symmetrical image about the real axis of those with the same value for ci. 
Thus, only circles for ci are shown in figure 6; the straight line represent- 
ing C£ = 1 is added to show the symmetry. 

The geometrical structure of the circles becomes obvious by considering 
the following properties (stated for the case of c* only): 

1. The center points of the circles given by a znil and to zm lie on the 
straight line passing through p^ and the origin. 1 


8 



2. The radius of the locus circles for ci = a and cj = 1/a is the 
same. The corresponding center points lie symmetrical to the point 
midway between p^ and the origin. The midpoint is on the locus 
circle for ci =1, which has degenerated to a straight line and 
intersects the real axis at 


a z = |p i | 2 /2o i (17) 

The line for ci = 1 is thus the perpendicular bisector of the line 
between p^ and the origin. 

3. The circles degenerate to a point at p^ or at the origin for 
values of ci going to either zero or infinity, respectively. 

The radii of the circles are given by the product between the 
distance of the center from the origin and the value of the 
corresponding magnitude. 

Considering the intersections of the ci and C 2 circles, each such 
intersection represents a value of |Cj_ z | , given by the product of 
ci and C 2 at that point. Vice versa, for a given point in the root 
plane, ci and C 2 can be determined and then | z | itself. 

Although the higher order loci curves of the magnitude of C^ z are given 
by the intersections of two families of circles, it still is rather tedious to 
derive those curves for each required case without great computational effort. 

A simpler interpretation can be obtained by encircling an area in the 
root plane for zeros that correspond to the required ranges of |c^ z |. That 
means the exact locus curve for the required value of | z | will be 

Figure 7 shows the upper half of 
the root plane (the lower half would 
show symmetrically the same), where 
for each of the three regions (I, II, 
III) separated by the locus curves 
for ci = 1 and C 2 = 1, an example of 
encircling is shown. Those examples 
are based on the following statements . 
(More details are given in appendix 
B.) 

Region I: (a) If 1 G-j^ z | is distinc- 

tively greater than 1, the circle of 
c 2 = |Cizl surrounds the area within 
which all conjugate complex zeros 
corresponding to that | C-j^ z | must 
lie (example, | C-^ z | > 2.0). (b) If 

I Ci z | is distinctively smaller than 
1, the circle of c i = | z | 


approximated by a surrounding region. 


i" 



Figure 7.- Examples of approximate 
zero locations for certain 
values of |Ci Z |. 
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surrounds the area within which all conjugate complex zeros corresponding to 
that | | must lie (example, \C± Z \ < 0.4). 

Region II: All conjugate complex zeros corresponding to a certain value 

of |c^ z | > 1 must lie within the union of the circles for c^ = |Cj[ z | and 
c 2 = l^izl "but outside the intersection of those circles (example, 

|Ciz| * 10). 


Region III: All conjugate complex zeros corresponding to a certain value 

of i c izl where 


i - <o i /|p 1 l) 2 < 



< i 


must lie within the union of the circles for c^ = | C-j_ z | and c| = |C^ Z | 
outside the intersection of those circles (example, | z | * 0.3). 


As the straight-line loci for cx = 1 and C 2 = 1 separating the different 
regions change with p-x , the circles encircling the zero regions are cut dif- 
ferently by those lines. However, the radii of the circles are proportional 
to | p-j_ | and do not change if pj_ remains constant in magnitude. 


For the locus of |Cj_ 2 | = 1, some other simplifying specifications can be 
made. In the relations given in appendix A, it becomes, evident that the locus 
for |C-x z | = 1 intersects the real axis at the same point that the loci for 
cx = 1 and c z - 1 intersect. The locus goes to infinity in its imaginary 
coordinate for 


a 

z 



2a. 

i 


(18) 


and at o = a . , u> is defined by 
z i 9 z 17 




(19) 


Notice that equation (19) can hold only if the locus intersects the real axis 
at the left of the pole p^. Otherwise, the relation for u z at a z = ( 1/2 )a ^ 
may be used, where 


1 2 1 

t - a = — 


+ ¥ a i 


2 

2 “i 


(20) 


Real zeros — Following the procedure worked out for the case of a pair of 
conjugate complex zeros (see fig. 6), it becomes evident that the intersection 
points of the circles for ci = const with the real axis represent locations 
for single zeros with | z | = c^. The same is true for those circles of 
C 2 = const, which leads to the same zeros and values of |C^ Z |. 
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For a single zero, | z | can be derived from the same set of curves used 
for the conjugate complex zeros , and no further set of curves for real zeros 
is needed. Real zeros thus are simply treated using a plot like figure 6 
separately for each. 



Figure 8.- Illustration of change 
in |Ci z | with variation of 
one real zero. 


How |C^ Z | changes as a conse- 
quence of shifting a simple zero 
along the real axis is further illus- 
trated by figure 8. As in figure 6, 
the zero location for |C iz | = 1 must 
lie at ttye point where the line for 
ci = 1 intersects the real axis (eq. 
(17)). From there |Ci z | increases 
toward infinity as the zero moves 
toward the origin. As the zero moves 
beyond the origin, J z J decreases 
monotonically to 1. 


As the zero moves in the opposite direction, the magnitude of Ci z 
decreases to a minimum value and approaches 1 again beyond the minimum point. 
The least achievable value of |Ci z | is dependent only on p^ , 


C. = 1 - 

12 min 



( 21 ) 


and the corresponding zero location is 


which is twice the value for the intersection point of the line of ci = 1 
(fig. 8). 


Interpretation of Pole Modifications 


The preceding discussion is concerned with C^ z , the zero-dependent part 
of the residues (eq. (6)), and the effect of zero changes on that part. The 
same concept can be used for the interpretation of Crp, 'the other part of the 
residues (eq. (6)), which changes only with modifications of the poles. Since 
the transfer function denominator is restricted to fourth order, the 
expression 


C. 

IP 


4 

n 


U=1 


_ -P-H_ 

- pu 
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can formally be discussed as that which would have been obtained as the recip- 
rocal of C^ z with three zeros, given the values of the poles p u . In the 
case of complex p.^ , the remaining three poles are used, in the same manner as 
before. However, since the "third zero" is the pole conjugate to p^ , its 
phase and magnitude are constant. The direct analog to the case of C-£ Z with 
a pair of zeros, from which the loci plots are formed, is C£p , the part of 
C-pp that does not contain the conjugate pole of p^. Thus, *the loci plots 
derived for C^ z can be -used for C*p as well, by relabeling the loci curves 
in a simple manner: 


and 



* -h 


CcJpJ 



XL 



( 22 ) 


( 23 ) 


APPROXIMATE MODAL EXPANSION METHOD 


Modal-Type Transfer Function Expansion (Two Zeros) 


The residue interpretation method discussed above is very appropriate for 
the interpretation of the effect of zeros on the time response, if the poles 
of the system under consideration are conjugate complex. In this case, with a 
fourth-order system, there are only two residues or time vectors that must be 
considered, each representing one of the eigenmodes . If instead of the conju- 
gate complex poles there are real ones , the number of independently moving 
time vectors increases (up to at most four). Pursuing the motion of these 
time vectors to determine their superposition at any time t is no longer 
simple enough. Even in the case of three independent time vectors the method 
loses much of its power. Therefore, in cases where some of the poles are 
real, a different method is proposed based on a modal-type transfer function 
expansion using typical pole locations obtained from longitudinal aircraft 
equations of motion. According to the transfer function given in equation (l) 
and its later normalization, the expansion is written for the case of a pair 
of zeros (b3 = 0) : 

r(J N(s) l+b lS +b 2 s 2 

D(s ) l+ais+a2S 2 +a3S 3 +a4S 4 

_ g I + 6 II + g III 

1+0t ll s+a 2l s2 1+a iII s+a 2II s2 (l+a 1;[ s+a 2I s2) (l+a 1I3; s+a 2II s2) 


- i e n £m 

DjCs) d ti (s) dST 


( 24 ) 
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The denominator of the third expansion term is that of the original 
transfer function * and the denominators of the two other terms are the charac- 
teristic polynomials of the eigenmodes of the aircraft called the phugoid (Bp) 
and the short-period motion (Bn)- These denominators will not be changed if 
one pair of the poles or even all poles become real. The t eigenmodes are con- 
sidered to be known in advance as their response characteristics. Thus the 
gains Bj and Bn represent a kind of modal expansion, and the Bin term 
should be incorporated in that expansion. With regard to the aircraft longi- 
tudinal motion, this can be done by an approximation (ref. 5). It will be 
assumed that Bm/D can be expressed using only the phugoid portion of the 
denominator. This is the case if (l) the values of Bi and Bin are no ^ very 
high and (2) the magnitude ratio between the resulting time vector of the sup- 
posedly aperiodic phugoid and that of the short period, both at time t = 0, 
is greater than about 10. Usually the magnitude ratio | p 3 —p 2 I / I P 1 — P 4 I i- s 
about 1, so the condition for the time vector ratio is met, if 


P3 1 > 3/pi * P2 


(25) 



Figure 9-~ Comparison between actual 
F(s) and second-order approximation. 


Figure 9 gives some evidence of the 
degree of coincidence that may be 
expected between the first and third 
expansion term time histories . 

With these assumptions , the 
expansion of the transfer function in 
equation (24) can be modified into 


F(s) = 


6 i * e m 


1 + a^s + a 2 i s " 
6 II 


1 + a „s + a s ‘ 

+ « 1];i a U 2II 


( 26 ) 


By this type of modal expansion, the superposition of only two modally charac- 
teristic time histories has to be made. The gain (6j + 8jjj) on the one hand 
and Sii on the other hand give the degree of dominance of one of the two 
portions. 


Zero Loci Due to the $ Values 

To assign zero locations to certain desired types of superposition in 
phugoid and short-period response given by , $n , and 8m , zero loci for 
constant 3 values can be plotted in the root plane. Equation (24) yields 
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N(s) = l+b 1 s+b 2 s 2 = 8 I (l+a 1II s+a 2II s 2 )+6 II (l+a 1I s+a 2I s 2 )+6 III 


( 27 ) 


so that Bj, Bjj, and Bju are given as quantities dependent on bi and b 2 , 
that is, on the zero locations, while the poles are considered as fixed: 
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These equations represent the loci curves. When written in terms of a z 
and to z for the case of conjugate complex pairs of zeros, they again are 
circles (fig. 10): 
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Figure 10.- 3j and 3 jj loci. 


with to 2 = a i + w z* case of pairs of real zeros with root plane axes for 
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they "become rectangular hyperbolas of the same algebraic structure: 
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With respect to 3p (or 3 jj) all hyperbolas intersect in one point, geomet- 
rically given by the intersection of the 45° line through the origin and the 
line o z = ai (or o 2 ). Figure 10 shows a plot of some loci curves for both 
conjugate complex (upper) and real (lower) pairs of zeros. The phugoid poles 
are fixed at pi = - 0.186 (l/sec), p 2 = -0.3 (l/sec), and the short-period 
poles at P 3,4 = -0.65 1 j 0.92 (l/sec). For simplicity, only loci for posi- 
tive values of 3j and 3 jj are shown. The 3m loci are also not plotted, 
because the value of 3m can easily be determined from 3q and 3 ji (eq. 
(28)). If the pair of zeros is located on either the phugoid or the short- 
period poles, 3m becomes zero. According to the approximation made in the 
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last section, the loci circles of constant 8m are almost coincident with 
those of 3 t , "but with different values assigned to them corresponding to 
equation (2o). In figure 10, the shading indicates regions of zero locations 
for which certain combinations in the signs of gj and gjj are valid. It 
appears, for example, that for 0 > cr z > ajj and gjj < 0, one of the zeros will 
he in the right half plane - a nonminimum phase effect. As a z approaches 
oil , the contribution of the phugoid to the adverse motion becomes larger and 
the phugoid portion is exclusively responsible for that kind of behavior when 
o z < ajj. 

As an additional interesting fact, there is no pair of conjugate complex 
zeros for 

(|) 2 TT 

oil 

^1 > " °j) ~ 

or 

u 2 ol 

3 II < " ^h^ii " V = 


For those values of either gj or $jj, the corresponding loci circles degen- 
erate into a point on the real axis at double aj or ai , respectively. On the 
other hand, the ranges in gj combined with certain values of gjj are also 
limited for certain regions of zero locations. For instance , in the example 
shown in figure 10, there is no gj > 0 together with gjj = 1.0, that will 
form a pair of conjugate complex zeros. 

For values of gj and gjj , as indicated in equations (32) and (33), the 
corresponding points in the a z ,u) z plane, each representing a pair of real 
zeros, lie on hyperbolas that do not reach the a z axis. As their conjugate 
axis is the o z axis itself, the point closest to the a z axis is given by 
the hyperbola vertex (fig. 10). With growing gj and decreasing gjj, the 
vertices of the corresponding hyperbolas are moving from the a z axis to a 
maximum displacement, where their vertices lie on the intersection with all 
other hyperbolas belonging to that loci family. The gj and gpj values for 
that case are 

8* = 2gj , g* x = 2g^ (34) 

It appears from figure 10 that there is only a small region of zero pair loca- 
tions in which significant contribution of the short-period modal response 
will be expected or equivalently where gqj is greater than 1.0. For these 
values , conjugate complex zeros must be located within the corresponding 
circle shown close to the origin. For pairs of negative real zeros the same 
is possible only in the narrow region between the 45°-line from the origin and 
the hyperbola extending the gjj = 1.0 circle in the lower half plane. 

Considering the case of single zeros in terms of 8j and. gjj, the defini- 
tion of g' and g* is also very useful. A single zero located at points 


(32) 

(33) 
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•where for pairs of zeros, Bp = gj and gjx = gjj (fig. 10), will have the same 
values for gj and gj j . This becomes evident, if the relations for gj, gjj, 
and gin are developed for the single-zero case from 

N(s) = 1 + bis = $ ] .(l + a 1IJ s + “ 2 u s2 ) + 3 i;[ (l + + a 2l s2 ) + 

(35) 

such that, with a = -1/bj , 

z 


6 I - 


°I 


■'ll 


III. 6 * 

a P II 
z 


(36) 


8 in - 1 - * 6 ii> 


These relations are very simple and do not require any plotting, as should be 
expected for the case of a single zero. 


Three Zeros 

The extension to the case of three zeros requires a slight augmentation 
of the transfer function expansion used in the last sections. An additional 
g coefficient must be generated to keep the main structure of the expansion 
and ensure the validity of the approximation that was made. The augmented 
expansion becomes 


F(s) = 


1 + b^s + b 2 s 2 + b 3 S 3 
1 + ajs + a 2 s 2 + a 3 S 3 + 8 , 4 s 1 


- . g II + SS II . g IH 

“ d^TsT d ti (s) D(s) 


Following the approximation made earlier, 


> 6 I + 6 III . 6 II + s6 XI 

F(S> • - DjCs )- * • - E r ~ ( ' s l 


(37) 


(38) 


This essentially means that, in contrast to the case of two zeros, only the 
short-period portion of the response will be modified by its derivative due to 
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the gain 8 jj. Neglecting the contribution of the 3 jj expansion term itself 
in the first place, 3 i, , and are all implicit functions of the third 

zero , and can be treated the same way in the case of a pair of zeros . The 3 
plots developed for the case of a pair of zeros (fig. 10 ) can be used for the 
interpretation of transfer functions with three zeros. In fact, 3j through 
3iu can be described by the relations of equation (28) by replacing (d qz 
and a z by (Dq Z and a z , as follows : 
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(40) 


which leads to the equation for the equivalent pair of zeros z { t 2 : 


where 


' 1,2 


= -A * /A* - B z 


A = 


-2z3a - 03 2 + co 2 T 

0 z oz ol 

2 (-Z 3 -2a + 2 a T ) 


B 2 = 


z 3 <d 


oz 


z 3 + 2a - 2 o t 
3 z I 


(4l) 


A root locus plot for z'1,2 as function of the parameter Q = Z3(a z -Oj)/2 can 
be derived (fig. 11). The root locus equation becomes: 


1 The third zero Z 3 is chosen as either the single real zero or arbitrar- 
ily from the three real zeros. 
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( 42 ) 


Thus , with increasing given mag- 
nitude of the third zero , the equiva- 
lent pair of zeros develops from a 
real pair (one at the origin) into 
the actual pair of zeros z{ £• 


Figure 11.- Equivalent pairs of zeros 
for use of Bp and Bn plot in 
the case of three zeros. 


EVALUATION OF OPTIMAL LONGITUDINAL CONTROL TIME RESPONSE 
WITH USE OF POLE-ZERO INTERPRETATION METHODS 


Although the methods described in the previous sections can be useful for 
analysis purposes in the frequency domain, they are used in this section for 
the evaluation of an optimal design. 

Linear optimal control theory optimizes the time response of a given sys- 
tem with respect to a quadratic performance criterion, weighting the states as 
well as the controls . In the case of the longitudinal motion of aircraft , the 
relation between the requirements on the time response and the weighting coef- 
ficients of the performance index is not strictly clear. Thus, only a hypo- 
thetical guess of the weighting coefficients can be made that has to be proved 
satisfactory by examining the actual performance of the corresponding time 
response and control laws. The dependence of the time response on the per- 
formance index may be studied effectively by the procedure outlined above. 

Here we study this dependence for the aircraft longitudinal control 
system. 


Problem Formulation 

The system investigated as an example is the regulator-controlled longi- 
tudinal motion of a conventional subsonic-type transport aircraft during land- 
ing approach (fig. 12). It can be described by the state equations: 

x = Fx + Gu = (F + GL)x (43‘) 

c 

where x is the state error vector 
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u c is the control vector 


Figure 12.- System "block diagram. 



and F, G, and L are matrices with constant elements. The F and G matrices 
are given as aircraft parameters , while the elements of L have to be deter- 
mined as the regulator feedback gains. This type of aircraft control is ade- 
quate when the flight path is prescribed by air-traffic control as a straight 
line, such that an equilibrium state can be assigned to it. Thus the task of 
the controller is to drive the errors of the states, caused by initial dis- 
placements x(0), to zero. Performing this task in the described setup, the 
regulator could consist of an automatic controller, or both automatic and 
manual in parallel. The total gains would be equal to L. If the gain por- 
tion due to manual control is kept inside the allowable range for pilot gains , 
there might be no further dynamical equalization required from the pilot to 
minimize the performance criterion chosen for the system. 


Performance Criterion 

To follow a given straight flight path corresponding to a certain equi- 
librium state, the task of the controller can be reduced to that of a regula- 
tor, mainly driving the magnitude of the error in the speed vector to zero. 


V c = commanded V due to 
desired path 

V a = actual V 



Figure 13.- Definition of the 
speed notations. 


The reduction of the magnitude of the 
speed error vector AV (fig. 13) is 
of. primary interest with respect to 
the overall behavior, irrespective of 
its direction. Thus, there is only 
one independent quantity describing 
the performance of the system output 
instead of the original four given by 
the components of the error state 
vector. This single error has to be 
weighted against the allowable effort 
in the controls with all controls 
having the same weighting. The mag- 
nitude of the error in the speed 
vector Vq can be expressed in terms 
of error state components , as shown 
in figure 13: 


1 For simplicity, the penalty of the error in height is ignored. Computa- 
tions have shown that this does not result in a big difference in the result- 
ing optimal control law. 
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| AV| 2 = u 2 + w 2 (44) 

With w = (u + |V C |) • tan y (45) 

and assuming small y as well as small u compared to I v qI > equation (44) 
becomes 

| AV | 2 = u 2 + ( | V c | • y) 2 (46) 

Thus , the minimum of the quadratic performance index 

I=f/" U(u 2 + |V C | V) + r(6 2 + « 2 + 6p) 3dt (1.7) 

defines an optimal regulator control law 

u = Lx (48) 

c 

with respect to the value chosen for the weighting ratio q/r, where q > 0 
and r > 0 (refs. 6, 7 ). 


Discussion of Resulting Optimal Transfer Functions 

Wow, suppose that L has been computed in the standard way. Talcing the 
Laplace transform of equation (43) with initial condition x(0), the response 
is given by a 4x4 matrix H, 

x(s) = [Is - (F + GL)] -1 x(0) = H • x(0) (49) 

whose elements can be considered as transfer functions with the error state 
components as outputs and the initial conditions as 6-function inputs. The 
elements of H are easily derived from equation (49): 

x^(s) det[ls - (F + GL)]^ 
h ij = Xj(0T = det [is - (F + GL)] 

where [.].. is equal to [.] in all elements except the ith column, which is 
replaced by the jth column of the 4x4 unity matrix. 

It will be shown that the main changes in the character of the time 
response take place in the range of q/r between 0 and 2 . 0 . Beyond that, the 
control effort becomes higher and higher, with no comparable significant 
changes in the time responses of the error states. For q/r = 2 . 0 , the 
phugoid poles are still conjugate complex, so that only the first described 
method is used. 
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The aircraft system has three controls and three states of particular 
interest , resulting in nine separate responses or transfer functions to be 
considered. Time histories of these states and controls for q/r values of 0 
and 2.0 are shown in figures ih through 17. The quantity y is also plotted 
to show the load factor. There are considerable changes in the time histories 
with increasing q/r. Figures 18 through 20 show the loci required for the 
application of this method. Figure 18 shows the poles and zeros for all the 
various transfer functions for q/r = 0 and 2.0. In addition, the phase loci 
are drawn for q/r = 2.0 for both p^ = Pi and p^ = P3. Figure 19 shows the 
poles and zeros , but also includes the magnitude loci associated with the 







Figure lU.- Transient responses due to 
initial conditions of original sys- 
tem (q/r = 0); (a) in 0, (b) in y, 
(c) in u. 


Figure 15 •- Transient responses of 
the state and control variables of 
the regulated system (q/r = 2.0) 
due to initial condition in 0. 
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Figure l6.- Transient responses of the state and control variables 
of the regulated system (q/r = 2.0) due to initial condition in y. 



t, sec 

Figure 17- - Transient responses of the state and control variables 
of the regulated system (q/r = 2.0) due to initial condition in u. 



Figure 18.- Phase loci for longitudinal response residues. 

phugoid residues again for q/r = 2.0. Figure 20 is similar to 19, with the 
magnitude loci drawn for the short-period residues . A number of the roots 
change only slightly with q/r, and some change very similarly for several 
transfer functions. 

In deriving the time vectors (or residues) by means of figures 18 through 
20, we use as an example the transfer function y/u(0) for q/r = 0, which has 
a second-order numerator and a fourth-order denominator.. As shown earlier, 
the time vector of the short-period mode Y3» for instance, will be determined 
from the residue C3 (at P3), which in turn consists of the product of Cjj 
and C3p (eq. (6)). The phase of C 3z ,<fi 3z is determined from P3 and zi as 
indicated in figure 21 (a), which shows the appropriate poles and zeros 
extracted from figure 18. The conjugate zero Z2 does not need to be included 
for this plot. The quantity |c 3Z | = ci • c 2 is obtained from figure 2l(b), 
which is abstracted from figure 20, using P3 and z 1 } z 2 . These quantities , 

^ 3Z and |c 3Z |, are then combined and plotted in figure 22 as C 3Z . A similar 
process, shown in figure 21(c) and (d), is used to derive psC^. 1 The phase 

1 Note that P3^ 3 p or P3C3 is used due to the 6-function input u(0) 
instead of a step input. 
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Figure 20.- Magnitude loci for longitudinal response short-period residues. 
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C 2 =8.4 



Figure 21.- (a) Phase locus and (b) 
magnitude loci for C 3Z of the 

transfer function y/u(0); q/r = 0. 



Figure 21.- (c) Phase locus and (d) 
magnitude loci for P3C 3 p of the 

transfer function y/u.( 0 ) ; q/r = 0 


of p3C 3 p, denoted as <P 3 p in figure 22, is obtained by using pi instead of 

z 2 (fig. 21 (c)). The angle of tt /2 must be added due to the contribution of 
p4 , the conjugate pole of P3, and the sign has to be changed (eq. (22)). The 
magnitude of P3C 3p is given by the product of | P3 | and |C 3 p| = 1/(01*02*03) 

which in turn is obtained from figure 21(d), using all poles pi through pi+. 
Again, phase and magnitude are combined to plot P3C3 in figure 22. 

The quantities C lz and pjC^ for the time vector Yi are obtained 
similarly using pi instead of P3. 



With P1C1 (P3C3) determined, 
the actual time vectors yi ( y 3 ) at 
t = 0 (fig. 22) are two times pjCi 
(P3C3). This same process is used tc 
determine all the time vectors shown 
in the subsequent figures . 

In the following, the time 
responses of the state components 0 , 
Y , and u due to initial conditions 
in the same states are discussed 
separately, referring to figures 18 
through 20. The time histories of 
the controls could be discussed in 
the same way. 

0/0(0) — The attitude is the main 
control variable in the longitudinal 
motion of aircraft. The requirement 
is to reduce the initial attitude 
error in a rapid, well-damped manner 
to zero. Without any control 
(q/r = 0), the time response behavior 
of the attitude due to initial dis- 
placement shows large overshoot and 
very little damping of the low- 
frequency mode, the phugoid (see 
fig. 14), whereas for optimal control with respect to the chosen performance 
index and increasing q/r the desired behavior can be achieved, as shown in 
figure 15. 

There are three real zeros but the overall tendencies in the change of 
the time response due to increasing q/r are caused mainly by the shift of 
the phugoid poles to higher damping and lower frequency. This shift will 
affect all time responses, because the poles are the same for all transfer 
functions for a given q/r. Therefore, for all transfer functions, there will 
be the tendency for Cj, the phugoid mode residue or half of the time vector 
01 at t = 0, to shift its phase into the imaginary axis (fig. 23), as the 
phugoid poles approach the real axis of the root plane. (A similar effect 
will occur for other transfer functions.) Because there is almost no change 
in the short-period poles with increasing q/r (also common for all transfer 
functions), the location of the zeros with respect to the phugoid poles deter- 
mines how far the phase shift in ©i will go and what change in magnitude of 
that residue will occur. 

As can be seen from figures 18 and 19, 1/Cp increases mainly in magni- 
tude and C 1Z in phase with higher values of q/r. The phase shift is 
chiefly due to the shift of the phugoid poles in relation to the zero that is 
closest to the origin. The change in the phase contribution of that zero 
amounts to almost ir/2 in the range of q/r from 0 through 2.0. Although 
there are also changes taking place due to the zero farthest from the origin, 
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Figure 22.- The time vector construc- 
tion from the residue quantities. 
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that zero has little effect on either 
the phase or magnitude of C 1Z . Thus, 
the tendency in C lz is simply, 
obtained by watching only the rela- 
tive location of the phugoid poles 
and the zero close to the origin. 

With respect to 63 or C 3 (the resi- 
due of the short-period mode), the 
change in C 3 p is due to the phugoid 
shift and the change in C 3Z is due 
to the shift of the zero far from the 
origin. These are opposing changes, 
such that the total change in O 3 is 
rather small. 


Figure 23.- Time vectors for the 

transient response of 0 / 0 ( 0 ) with 
q/r = 0 and 2 . 0 . 


The conclusion is that the 
effect of the phugoid mode is con- 
siderably reduced by the large phase 
shift in 0 i toward the imaginary 

axis and the high phugoid damping. On the other hand, choosing a well-damped 
location for the phugoid and fixing one of the zeros , a given range of the 
Ci z -phase loci can be assigned to the location of the remaining zeros to 
achieve the desired behavior. 


y /Q(0) and u/dCO ) — Although suppression (decoupling) of the responses in 
the flight-path angle y and the speed u due to initial displacement in the 
attitude might be desired, a certain amount of temporary displacement in both 
y and u is required to avoid excessive costs in the controls. That is, a 
very effective lift control would be needed, with simultaneous speed control, 
to compensate for the initial lift change due to 0 ( 0 ). 

As can be seen from figures l4 and 15, at least some improvement is made, 
particularly in the speed response, by minimizing the previously described 
performance index. Again, the main reason is the shift of the phugoid poles 
toward the real axis of the root plane, which leads to high damping and fixes 
the phase of yi and close to the imaginary axis for both the y and u 

response (fig. 24). 

To determine the actual phase of yi and Y 3 j one can make use of the 
results obtained by examining the zeros for 0 / 0 ( 0 ) that are located in close 
vicinity to both existing zeros of y/0(O). Therefore, the phase difference 
between the residues of the responses in attitude and flight-path angle is 
caused by the phase effect of the medium zero of 0/0(0). Because this zero 
does not change with q/r, the phase difference between the residues of 
0 / 0 ( 0 ) and y/ 0 (O ) is also independent of q/r. 

Response due to y(0) ^ 0 — As shown in figure 17 , the response to an ini- 
tial displacement in the flight-path angle is good for the case of q/r = 2 . 0 . 
The flight-path angle itself settles down to small values very soon without 
evident oscillations , and the coupled motion in pitch and speed is almost neg- 
ligible. Also, the load factor due to the rate of change of the flight-path 
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Figure 24.- Time vectors for transient 
responses of y/0(O) and u/0(o) 
with q/r = 0 and 2.0. 


Figure 25 •- Time vectors for trans- 
ient responses of 0/y(O), y/y(0), 
and u/y(0) with q/r = 0 and 2.0. 


angle does not appear to yield too high values . The plot of the time vectors 
of the response in 0 is shown as part of figure 25* The short-period vector 
is phased with respect to the phugoid vector such that. while rotating and 
before decaying, their real parts almost cancel, resulting in low coupling 
between 0 and y(0). This configuration of the vectors thus shows a highly 
decoupled transfer function, and a similar configuration would be desirable 
for other transfer functions. Comparing this configuration with that of 
y/0(O) (fig. 24), for instance, a smaller phase shift for Y3 would produce 
less coupling. This could be obtained by locating the zeros on the corre- 
sponding phase locus curve. These phase values are easily obtained from fig- 
ure 25 since only one zero is involved in the 0/y(O) response. 

The pair of conjugate complex zeros of y/y(0) is almost coincident with 
the zeros of u/y(o) and y/u (0). Thus, the difference in the behavior of 
these later two in comparison to y/y( 0) is given by the real zero of y/y( 0 ), 
which again is almost independent of q/r. Observing the magnitude ratio 
• between the residues Yi and Y3 in figure 25, one can see that increasing q/r 
results in a high gain for yj relative to Y3, which gives a very slow long- 
term dynamic behavior. 

Response due to u(0) ^ 0 — According to figure IT, the direct response in 
the speed again appears satisfactory for moderate values of q/r. There is no 
fast response required for the speed so long as the initial error returns to 
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Figure 26.- Time vectors for transient 
responses of 0/£i(O), y/u(0), and 
u/u(0) with q/r = 0 and 2.0. 


q/r 


zero in a monotonic "but not too slug- 
gish manner. This is achieved by the 
high damped, practically pure phugoid 
motion (fig. 26) due to the pair of 
conjugate complex zeros very close to 
the short-period poles (figs. 18-20). 
On the other hand, the coupling of 
the velocity to flight-path angle and 
the pitch angle is fairly large, giv- 
ing somewhat unsatisfactory responses. 
Although there is a considerable 
improvement in the overall dynamic 
behavior for q/r = 2.0, compared to 
the case without regulator control 
(q/r = 0), the gains K of the 
transfer functions, especially of 
y/u(0) , are still high. 


The effect of the zero of 
0/u(O), which has a large decrease in 
value as q/r varies from 0 to 2.0 
(figs. 18-20 ), is to cause a large 
phase shift of the short-period time 
vector (residue) (fig. 26). Also, 
the phugoid vector moves across the 
imaginary axis before merging with it 
for q/r > 2. Considering that both 
the phase shift of the residues and 
the high damping of the phugoid are 
responsible for the improvement in 
less than 2.0 might do more to minimize the peak in 


the 0 response, a 

the 0-time history. It can be deduced from figure 26, by considering the 
motion of the time vector 0i for values of t greater than 0, that a smaller 
q/r will reduce slightly the peak in 0, but increase the long-period 
response. This result is shown on figure IT- Note that the other responses 
will also change, which must be considered. 


Abstract of the results — The simple performance index, weighting only the 
error state in the flight-path angle and the airspeed against the controls , 
using one parameter for each, provides in the main satisfactory results. For 
q/r ~ 2 most of the transient responses are satisfactory, but two of them are 
only marginally acceptable. 


As shown in figure 17, improvement in those time responses would require 
a different performance index. Since the error in y is relatively more 
annoying than that in u, a lower weighting of u compared to y would be 
desirable. Although this improvement is not carried out in this study, it 
would require at least one additional weighting parameter, and two parameter 
patterns of root loci drawn coincident with the optimal results for different 
levels of weighting in both parameters. The methods presented here could 
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still be used to aid in the determination of the elements of q and r that 
would give the desired time responses. 


CONCLUDING REMARKS 


Two methods are described that provide geometrical insight into the 
effect of the zeros of transfer functions as well as poles, on the correspond- 
ing time histories. The main advantage of these methods lies in their simplic- 
ity compared to other available procedures. They differ in their limitations, 
which are essentially those on the possible regions of the system pole 
locations . 

The first method, which uses the geometrical interpretation of the resi- 
dues, loses much of its power if the poles become real, while the other 
method is restricted to certain separations of the poles, regardless of 
whether they are real or conjugate complex. In almost all cases, at least one 
of the methods can be used effectively. This is especially true when applied 
to aircraft transfer functions, as shown in the discussion of optimal regula- 
tor control for the longitudinal motion. 

Both methods are developed for the type of transfer function with four 
poles and two zeros along with an extension to three zeros. Further extension 
to allow for additional poles and zeros is possible, but will be very 
complicated. 

The main use of the methods should be for analysis purposes. To inter- 
pret the effect of zero changes, the poles are assumed fixed, and changes in 
the poles are discussed with the zeros unchanged. Thus, only if something is 
known about the sensitivities of the pole locations due to zero changes and 
vice versa (ref. 8) can the methods be helpful in the synthesis process. 

The methods are applied to a study of the longitudinal motion regulator 
of an aircraft. The performance of this regulator system is studied using a 
particular performance criterion with minor variations. This criterion mini- 
mizes the state errors effectively using only reasonable control effort. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 9^035, Feb. 3, 1972 
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APPENDIX A 


EQUATION FOR THE MAGNITUDE LOCI OF C iz (TWO ZEROS) 


The exact equation for the magnitude loci of Ci z is simply derived from 
equation (l4) hy substituting in the expressions for ci and C 2 the values of 
the conjugate complex zeros, zi,2 = °z * j u z : 


C 

1Z 


= Cl • c 2 


[(a - a.) 2 + (to - to. ) 2 ][(a - cr.) 2 + (to + to.) 2 ] 

zi zi zi zi 


(a 2 + to 2 ) 
z z 


(Al) 


The coordinates of the pair of conjugate zeros generate fourth-order terms, 
such that the loci of iC-r^l in the root plane "become higher order curves. 
For the special case of \C±z\ = 1, which might he of particular practical 
interest, equation (Al) reduces to 


4a 3 a. 
z 1 


2a 2 (oo? + 3a?) + ka a. 
z 1 1 z 1 


p. I 2 + ka co 2 a. 

1 1 Z Z 1 




0 


(A2) 


From this equation, the simple formulas (lT) through (20) are derived for 

specific points on that curve hy substituting the specific values for a z or 

a) . 
z 


For a pair of real zeros, equation (A2) changes to 

Ua 3 a. - 2 a 2 (to? + 3a?) + ka a. |p. | 2 - ka co 2 a. + 2co 2 (a? - co?) 
zi z 1 i z i lir i 1 zzi zi 1 



(A3) 
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APPENDIX B 


ENCIRCLING OF ZERO REGIONS FOR GIVEN VALUES OF |C. 

1 1Z 


Practicable approximations can be made to encircle zero regions for given 
values of | C-j^ z | rather than to define the exact loci curves. This appendix 
gives some further details on those approximations. 

According to equations (15) and (l6) for the loci circles of Ci and C 2 , 
figure 27 illustrates the relation between the location of the center points 
of those circles and the corresponding values of ci and C£. It shows that on 
the string of center points of the ci circles, the corresponding value of ci 
is 

c? ■ 1 + - > 1 (Bl) 

1 n 

if the center point lies at a distance n|p^| downward from the origin and 

c? = 1 - - < 1 (B2) 

1 n 


if the center point lies n | pi | upward from the origin. The centers of the 
C 2 circles are symmetrical about the real axis with the corresponding ci 
centers . 



Figure 27.- Relation between values of 
c i ,C 2 and the location of the loci 
center points . 


In looking for the intersections 
between ci and C 2 circles, it is 
evident that not all ci circles 
intersect a particular C 2 circle and 
vice versa. Therefore, considering 
only the upper half of the root plane 
for symmetry reasons , there are three 
regions divided by the straight-line 
loci of ci = 1 and C 2 = 1 (fig. 7), 
such that (l) region I includes the 
intersections possible between the 
circles of ci < 1 and C 2 > 1, (2) 
region II includes the intersections 
possible between the circles of 
Ci > 1 and c 2 > 1, and (3) region III 
includes the intersections possible 
between the circles of ci < 1 and 
c 2 < 1. Encirclement of zero loca- 
tion areas due to certain values of 
|Ci z | is given below for each 
region. 


3b 



Region I 


If one considers a possible intersection between circles of cj < 1. and 
c 2 > 1 (n > 1, m > 0) , the corresponding magnitude of C iz can be written 
with use of equation (lU): 

i°i«l 2 -( 1 -s )( 1 + ii) (B3) 

or as an explicit relation for m as a function of n with the parameter 
|c. | , as illustrated in figure 28, 

X z 


m = 


n - 1 


(|C iz | 2 - 1) • n + 1 


(Bit) 



The figure shows that for values of 
|C^ Z | distinctly smaller (larger) 
than 1, n (m) does not go over a cer- 
tain small value. Therefore, the' 
range of circles in ci (c 2 ) corre- 
sponding to those (Ci z | values is 
also very limited. Thus, for a cer- 
tain value of |Cj_ z |, the correspond- 
ing zeros and all those for smaller 


'1Z 1 


must lie 


(greater) values of 
within the largest possible circle 
for ci (c 2 ), which is that for 


C 1 = |Ci z |?c 2 = |c iz |) (see fig. 7). 
The range of possible combinations of 
m and. n is not only limited as shown 
in figure 28; for higher values of n 
and m, the circles might not inter- 
sect such that no locus point for 
|Cizl can be defined in the root 

plane. The upper limit for m and n can be easily described by the following 
condition, which depends on the location of the pole pj_: 


I C 1 


c 2 = 2 


Oi. 

1 


I Pi 


(B5) 


Figure 28 itself is valid for all locations of p^ . 


Region II 

Within this region the relation between the location of the circles for 
ci > 1 and c 2 >1 (n>0,m>0) and | C-j_ z | itself is 
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(B 6) 


or 




in 


n + 1 

(| c . z | 2 - l)n - 1 


(B7) 



Figure 29.- |Cj. 2 | as a function 
of m and n (region II).. 



Figure 30.- Approximate region for 
I Ci Z | locus (region II ). 


as illustrated in figure 29 , which 
again is independent of the location 
of the pole p^. The curves are sym- 
metrical with respect to the straight 
line m = n. Calling n = n* and 
m = m* for the point m = n , for each 
curve it can be stated that for all 
n > n* (lower half of the curve) 
m < m#, and also for all m > m* 
(upper half of the curve) n < n*. 
Therefore, each point on the lower 
half of the curve represents a zero 
location that must be inside the C£ 
circle given by m = n* and outside 
the ci circle given by n = n*. A 
corresponding remark can be given for 
the upper part of the curve . Thus , 
the loci for | z | can be located 
in the root plane as: (l) being 

within the union but outside the 
intersection of the circles for ci 
and C 2 corresponding to n = n* and 
m = m*, respectively, where 



and (2) intersecting the real axis at 
the same points, where the circles of 
(l) are intersecting (fig. 30). 

From equation (B6), it follows 

that 



As in region I, there is an upper 
limit for m and n , dependent on p£ , 
which is defined as follows: 


36 





U). 

1 


Cl - c 2 =2 


I Pi 


(B9) 


This condition can he used, of course, for further encircling. 


Region III 

Analogous to the approximation applied to region II, for region III, with 
ci < 1 and C 2 < 1 (n > 1, m > l) , 


or 



m 


1 - n 

(|C. z | 2 - l)n + 1 


(BIO) 


(Bll) 


which is illustrated in figure 31. 
Again, the loci for |Cj. z | in the 
root plane can be located as within 
the union but outside the intersec- 
tion of the circles of ci and C 2 
corresponding to n = n* and m = m* , 
respectively, with 


n 


Figure 31.- |Ci z | as a function of 
m and n (region III). 


C 1 - c 2 - J ^iz 


and 



(B12) 


(B13) 



The upper limits on n and m in figure 31, dependent on the location of , 
are given by 


to. 

| c ! + c 2 | =2 — — (BlU) 

IPil 
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